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Abstract 

A wide class of materials acquires magnetic properties due to particle interactions 
through exchange forces. These can be atoms and molecules composing the system 
itself, as in the case of numerous magnetic substances. Or these could be different 
defects, as in the case of graphene, graphite, carbon nanotubes, and related materials. 
The theory is suggested describing fast magnetization reversal in magnetic systems, 
whose magnetism is caused by exchange interactions. The effect is based on the cou¬ 
pling of a magnetic sample with an electric circuit producing a feedback magnetic field. 
This method can find various applications in spintronics. The magnetization reversal 
can be self-organized, producing spin superradiance. A part of radiation is absorbed 
by a resonator magnetic coil. But an essential part of radiation can also be emitted 
through the coil sides. 
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1 Introduction 


Magnetic materials, whose magnetism is cansed by exchange interactions, form a very wide 
class of magnets of different sizes, inclnding the samples of nanosizes, which hnd a variety 
of applications |T] . A novel class of exchange-interaction magnetic materials is the graphene 
family with defects, inclnding graphene flakes and ribbons, graphite, and carbon nanotnbes 
HU. Magnetic materials hnd nnmerons applications in qnantnm electronics, for example, 
in spintronics, information processing etc. 

In the present paper, we concentrate on two interconnected featnres of magnetic mate¬ 
rials: (i) First, we study the way of fast magnetization reversal in these materials that is 
a property crucially important for spintronics and information processing. We show that 
by coupling a magnetic sample to a resonant electric circuit it is possible to achieve fast 
magnetization reversal. The details of the reversal can be easily regulated by varying the 
system parameters, (ii) Second, fast magnetization reversal should produce magneto-dipole 
spin radiation in radio-frequency or microwave region. However, a large part of this radiation 
would be absorbed by the coil of the coupled electric circuit. We aim at studying whether 
some part of the maser radiation could be emitted through the coil sides. The possibility of 
this effect would essentially widen the region of applicability of such magnetic materials. 

The fact that coupling a resonant electric circuit to a spin system could essentially in¬ 
fluence spin dynamics is the essence of the Purcell effect [5]. A detailed theory of spin 
dynamics, employing the Purcell effect, has been developed for nuclear spins [6HIQ], mag¬ 
netic nanomolecules [nHl3], and magnetic nanoclusters [HHIZ] (see also the review arti¬ 
cle m)- The characteristic feature of all these materials is that their particle interactions 
are described by dipolar forces. Also, only the total radiation intensity has been consid¬ 
ered. However, in the typical setup, the magnetic sample is inserted into a coil of a resonant 
magnetic circuit, so that an essential part of radiation is absorbed by this coil. In order to 
understand whether some part of magneto-dipole radiation could be emitted through the 
coil sides, it is necessary to study the spatial distribution of the radiation. 

The main novelty of the present paper is twofold: (i) We consider the class of magnetic 
materials with exchange interactions. This class is widespread and important, including a 
number of known magnetic materials as well as new nanomaterials, such as graphene with 
defects, (ii) We study the spatial distribution of radiation in order to conclude whether a 
sample inside a coil could serve as an emitter of radiation passing through the coil sides. 


2 Magnets with exchange interactions 

The Hamiltonian of a system of N particles is the sum 

N 

H = ■ ( 1 ) 

i=i 

The first term is the Hamiltonian of an anisotropic Heisenberg model, 

lJ‘j {Sts’ + S!S>) + i.jS!S]] 


( 2 ) 


describing particles with exchange interactions. The second is a Zeeman term, with po being 
magnetic moment, and with the total magnetic held 
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(3) 


consisting of an external magnetic held Bq along the z axis and a feedback held H, along 
the X axis, caused by the magnetic coil of an electric circuit. The sample is inserted into the 
coil, with its axis along the x axis. The feedback held is dehned by the Kirchhoh equation 
that can be written [8l|9l|18] as 

— —h 2'yH + u^ H{t ) dt = —dvr —;— . (4) 

dt J Q dt 

Here 7 is the circuit attenuation, oj is the circuit natural frequency, and the ehective elec¬ 
tromotive force is due to the moving magnetization of the sample. 


N 

^ E (S’) . (5) 

i=i 

where I 4 is the coil volume and angle brackets imply statistical averaging. 

The 2 ; axis is assumed to be the axis of the easy magnetization, so that lij should be 
larger than Jij. If exchange interactions are due to electrons, then /tq < 0. Therefore, if 
Bq > 0, then the equilibrium spin value So of a particle is negative, Sq < 0. In this case, the 
Zeeman frequency is positive, 

1^0 = —ho-Bo > 0 . (6) 

The ehective particle spin S can be arbitrary. The ladder spin operators 

Sf = SJ ± iS] 

satisfy the commutation relations 

[s+, 7] = 2%7. [ 7 , 7] = 2%7 . 

In terms of the ladder spin operators, the exchange Hamiltonian (2) is 

5 E SiSts- + I„S’S‘) . (7) 


Writing down the Heisenberg equations of motion, we compliment them by the atten¬ 
uation 7 i, caused by spin-lattice interactions and 72 , due to other spin interactions, e.g., 
dipole interactions that usually are smaller than the exchange interactions. I that way, the 
equations of motion read as 


± 


dt 


Sf. H 
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dt 


S], H 
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( 8 ) 


Accomplishing the related commutations, we get the equations for the transverse spin. 


dt 


t Y, (hiSrS‘ - J„S‘S-) - !u)„ 7 - Vo/fSr - ^,_S- , 


(9) 


and for the longitudinal spin 

d ? ? 

= 2 ^ ~ + 2 (St - Sr) - 71 (Sf - So) . (10) 











3 Stochastic mean-field approximation 


We shall be interested in the evolution of the averaged quantities characterizing the transition 
function 


u = 


1 

Ws 
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E . 
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the coherence intensity 


and the spin polarization 
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Averaging equations (9) and (10), we invoke stochastic mean-field approximation [19], 
replacing spin pair correlators as 


{SfS^) ^ {St){S^) + {Sf)5S^ + {S^)5S: , 


(14) 


where i j and 5Sf is treated as a stochastic variable, such that its stochastic average be 


zero: 


{{iS“))=0 


(15) 


Then the averages {Sf) become functions of the stochastic variables. 

Also, we introduce the composite stochastic variables 

s ^ - h,ssp , { s (k,SS- - Ji,SS-) 

j{N) 

- isj) . 

j{N) 

And we dehne the effective anisotropy 


(16) 
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We assume that {Sf) does not depend on the index enumerating spins and that the variables 
(16) also are not index-dependent. 

In this way, from equation (9), we obtain the equation for the transition function. 


dzL 

— = -i{uo + ^0 - SAJz - i'y 2 )u - fiiaoH - ^)z , 


and for the coherence intensity, 
dw 


dt 


= -2^2W + i {noH - C) ZU - i (/xo-ff - 0 


(18) 


(19) 



While equation (10) yields the equation for the spin polarization, 

(17 1 1 

- = ( 20 ) 

where ( = Sq/S. The initial conditions 

Mo = ^i(O) , Wo = m;(0) , zo = z(0) , 
compliment the above equations. 


4 Resonator feedback field 

The Kirchhoff equation (4) can be represented |Hlll3] as the integral equation 

H = —47r [ G{t — dt' , (21) 

Jo 


with the transfer function 


r _ y _ 

G{t) = cos{ut) — — sin(a;t) 
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and where 


lioNS d 

rrin-^ = —7TZ - - [u + u) . 
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The longitudinal and transverse attenuations are assumed to be small, as compared to the 
Zeeman frequency, 

(23) 


— < 1 , — < 1 . 
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The coupling of the magnet with the resonant electric circuit leads to the appearance of the 
coupling attenuation 

N 

7c = y ■ (24) 

The latter, together with the circuit attenuation, are small, as compared to the resonator 
natural frequency, 

(25) 
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The stochastic variables (16), because of condition (15), can also be treated as effectively 

small. 

To be efficient, the resonator has to be tuned to the Zeeman frequency, so that the 
resonance condition 

|A| 

—^ Cl (A = (u — (Uo) (26) 


00 


be valid. Also, the magnetic anisotropy has to be small, with the anisotropy parameter 

A = -< 1 . 

OOo 


( 27 ) 


Otherwise the sample magnetization would be frozen. 







In this way, the feedback equation (21) can be solved iteratively, taking for the initial 
approximation u ~ Mq exp(—with 

a ;5 = Wo — SAJz = Wo(l — Az) . (28) 

Then in the hrst iteration, we get 

= i {uip — u*ip*) , (29) 


with the coupling function 


1 — exp{—i(w — us)t — jt} ^ 1 — exp{—i(w + us)t — yt} 


7 + i(w - Us) 


7 - i(w + Us) 


(30) 


In view of the above conditions, the hrst, resonant, term of the coupling function prevails, 
so that 

1 - exp{-iAst --ft} 

= -fcOJs - 


7 + iAs 

with the effective dynamic detuning 

As = u — Us = A + uqAz . 
Dehning the dimensionless coupling parameter 


9 = 


7c^o 

772 


(31) 


(32) 


(33) 


the real and imaginary parts of the coupling function can be written as 


R*’/’ = 9 .7,72 (1 - 212) p 


7 " + A| 
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and, respectively. 


ImV' = -g (1 - -47) P 
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cos(A5't) — — sin(A5t) 
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cos{Ast) + sm{Ast) 
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(34) 
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(35) 


Note that in the case of resonance, when A = 0, and under weak anisotropy A 1, 
we have As —)■ 0. Then the imaginary part (35) is close to zero, and the real part can be 
simplihed to 

ReV* ( 772(1 — Az) (1 — e”'*'*) . 

The latter form can be employed only when the anisotropy is not strong, so that A 1. 


5 Scale separation approach 

Substituting equality (28) into the evolution equations (18), (19), and (20), we have the 
equations for the transverse function, 

du 

— = -i{us + ^o)u - (72 - zil))u - zil)*u* + i^z , (36) 



the coherence intensity, 


dw 

dt 


= — 2(72 — az)w + i {u*^ — ^*u) z — z , 


(37) 


and for the spin polarization, 


^ = -aw - ^ {u*^ - Lp*u) + ^ [i}u^ + ip*{u*f] - 7 i(z - C) • 


Here the notation 


a = - (%()* +%()) = Rei/j 


(38) 


(39) 


is introduced. 

We solve the system of equations (36), (37), and (38) by resorting to the scale separation 
approach p,l^[^IT3 lfT8t lT9] . The function u is classihed as fast, while w and z, as slow. First, 
we solve equation (36) for the fast variable, keeping the slow variables as quasi-integrals of 
motion and taking account of the resonance condition (26). This gives 

u = uq exp < -{icos + 72 - - i / dt' > + 


(40) 


+ y ^{t')exp + ^2- zijj){t-t') -i J ^{t") dt" j dt' . 

Stochastic variables, in view of condition (15), are zero-centered, 

mm = mm = mt))) = 0 ■ 

The stochastic variable ^0 is real, while ^ and (j) are complex-valued. Therefore 

mmm = mitMm = {mm))) = o. 

For the non-zero pair stochastic correlators, we set 

aemm = - f) , = 2j,s(t - f), 


(41) 


(42) 


(43) 


with 73 playing the role of an attenuation caused by stochastic fluctuations. 

Substituting expression (40) into equations (37) and (38) and averaging the latter over 
fast temporal oscillations and stochastic variables, we come to the equations for the guiding 
centers, describing the coherence intensity. 


dw 

dt 


- 2(72 - az)w + 2'ym 


and the spin polarization. 


dz 

dt 


—aw — 73 ^ — ji{z 


Recall that a is given by equations (34) and (39). 


0 - 


(44) 


(45) 


6 Triggering spin waves 

Stochastic variables are the fluctuations that trigger spin motion. To clarify the nature 
of stochastic variables, let us show that these correspond to spin waves describing spin 
fluctuations around the average local spin polarization = z{t). It is worth recalling that 
spin waves can be well defined for nonequilibrium systems |20l[2T]. 

Spin fluctuations are defined as small oscillations around the average spin, which can be 
represented as 

Sf = (Sf) + 6S^ . (46) 

Substituting this into the equations of motion (9) and (10), we omit, for simplicity, the 
attenuations, keeping in mind the initial stage of the process, when time is yet much shorter 
than the relaxation times. Separating the zero-order equations, we get 

^(S-) = -..s(Sr). 4(y>=0. (47) 

And to first order, we have 

4 ss- = + (Spio - (spo . 

4 i ((S-) 4 =- - (Spp . (48) 

Keeping in mind that z is a slow variable, from equations (47), we And 

(S-) = , {S]) = zS . (49) 

In order to stress the role of fluctuations, let us set Uq = 0. Then 

5S- = S- , 5S; = 0 {uo = 0) , (50) 

where we take into account condition (15). The first of equations (48) yields 

^ S~ = -iuoS~ + izSi . (51) 

Let us employ the Fourier transformation for the ladder spin operators, 

■S/ = E , S,- = 15 ; S-er<^- , 

k j 

and for the exchange interactions, 

Jij = , {rp = n - rj) 

with the similar transformation for Ip. For the stochastic variable we get 


Then equation (51) reduces to 
with the spin-wave spectrum 


d_ 

dt 


c— • c- 
— , 


— ^0 + zS{Jk — Iq) ■ 
For long waves, when A; —)■ 0, we obtain 










This demonstrates that stochastic fluctuations are nothing but spin waves. 
For what follows, it is useful to keep in mind the restriction 


w + z = 


N 


NS 




i=i 


< 1 


(52) 

(53) 

(54) 


(55) 


which is necessary to take into account when setting initial conditions for equations (44) and 
(45). 


7 Qualitative classification of regimes 

Suppose that the initial state of the magnet is nonequilibrium, such that the external mag¬ 
netic held is directed along the ;7-axis, together with spins. Spin waves trigger spin motion, 
forcing spins to move, which generates the resonator feedback held acting back on spins and 
collectivizing their motion. The spins tend to reach an equilibrium state, reversing from 
positive to negative values. The overall dynamics can be classified onto several qualitatively 
different stages. Below, we give this qualitative classification, setting for simplicity zero 
detuning A = 0 and considering the limit of negligibly weak anisotropy A —)■ 0. Thus, we 
set 

A^O, A^O , As^O. (56) 

Then the coupling function (39) reduces to the simple form 

a fl'72 (l - • (57) 


7.1 Chaotic stage 

At the beginning, before coherence in spin motion sets up at time tcoh, the process is yet 
chaotic. At short times in the interval 

0 < t < tcoh , (58) 

the coupling function is yet very small, such that 

a 72 , a 73 . (59) 

Then the linear in time solutions of equations (44) and (45) are 

W Z^Wo + 2 ['JsZq - 72Wo) t , 


z ^ Zo- [(2:0 - C) 7 i + 73^0] t ■ 


(60) 



This regime lasts till the time, when the coupling function grows, so that 


(T-Z 72 {t tcoh) ) 

(61) 

which dehnes the coherence time 


tcoh 1^(1 • 

7 9Zq-1 

(62) 

At large coupling parameter, this gives 


tcoh = {. 9 Z 0 > 1 ) • 

IcWoZo 

(63) 

7.2 Coherent stage 


After the coherence time, spins become well correlated by means 
held. In the time interval 

of the resonator feedback 

tcoh < t < T 2 , 

(64) 

where T 2 = I/ 72 , the coupling function is large. 


CK > 7i , a > 73 , 

(65) 

reaching the value 


a ~ 5f72 (72 < 7 ) • 

( 66 ) 


The spin-lattice attenuation is usually small, such that 71 72 and 71 73 . Then 

equations (44) and (45) read as 


dw dz 

— = - 272(1 - gz)w , — = -g-i 2 W . 


The exact solutions of these equations give the coherence intensity 
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sech^ 


\9i2J 


t-to 


and the spin polarization 


2 : =-^ tanh 

912 




t-to \ 

/ 9 ' 


(67) 


( 68 ) 


(69) 


Here to and Tp are the integration constants that are dehned by the values w(tcoh) = Wcoh 
and z(tcoh) = Zcoh- This gives the delay time 


to = tcoh + y 


Ip + lg 
Ip-lg 


and the pulse time 


Tp = 


1 

Ip 


where 


H = 


1 + Wl + 4 


^72 \ 

I Wcoh 

% / 


7p = l2(gZa,h - 1) . 


(70) 


(71) 


(72) 



Under a large coupling parameter and small coherence intensity, at the coherence time tcoh, 
we get 

Ig - gi2Zcoh {gZcoh > 1 ) , 


Ip ^ gi2 \J + Wcoh {Wcoh -C Zcoh) ■ 
The maximum of the coherence intensity occurs at the time to, being 


(73) 




Wcoh + 




At this moment of time, the spin polarization is 

^(to) = - • 

g 


(74) 


(75) 


7.3 Relaxation stage 

After the transverse decoherence time T 2 , there is the relaxation stage in the interval 

T 2 < t < Ti , (76) 

where Ti = I/ 71 . For t S> to, fhe coherence intensity relaxes as 


w ~ 4tc(to) exp 

and the spin polarization tends to 





-Zcoh H-1-2 

g 





(77) 


(78) 


7.4 Quasi-stationary stage 

In the limit of very long times, when 


Ti < t < cx) , (79) 

the solutions exhibit small oscillations around the fixed points w* and z* given by the equa¬ 
tions 

72(1 - - 73 (^*)^ = 0 , 5 - 72 ^*+ 732 ;*+ 7 i(^* -C) = 0 . (80) 

7.5 Punctuated superradiance 

The coherent stage of radiation corresponds to the regime of superradiance, which can be 
pure, if Wq = 0, or triggered, if Wq > 0. It is also possible to realize the regime of punctuated 
superradiance [18], when in the process of spin dynamics either the external magnetic held 
is reversed or the magnetic sample is rotated, so that to reproduce the initial nonequilibrium 
conditions. In this case, radiation exhibits a series of coherent pulses, with the temporal 
intervals that can be regulated. 



8 Numerical investigation of dynamics 


We solve equations (44) and (45) numerically, with the coupling function dehned in equations 
(34) and (39). The case of exact resonance is assumed, with oj = uq. Then A = 0 and 
As = uAz. The chosen initial conditions correspond to a purely self-organized process, in 
the absence of imposed initial coherence, so that Wq = 0 , and when the initial polarization 
Zq = 1 dehnes a strongly nonequilibrium state, since the equilibrium polarization of a single 
spin, under the given setup, corresponds to C = — 1 - 

In the presence of a resonator feedback held, spin reversal happens much faster than 
the homogeneous transverse relaxation time T 2 = 1 / 72 - We measure all attenuations and 
frequencies in units of 72 . The spin-lattice attenuation 71 is usually much smaller than 72 . 
Taking this into account, we set 71 = 0.001. The dynamic attenuation 73 is assumed to be 
of order of 72 , which, in units of the latter, implies 73 = 1. Time is measured in units of T 2 . 

First, we consider the more general form corresponding to expression (34), which reads 


as 


_ 7 ^ 72(1 - Az) 

'y^ + iooAzy 



cos{uAzt) 


u 

7 


Azsm.{ujAzt) 



The temporal behaviour of the coherence intensity and spin polarization, for different system 
parameters, is shown in Figures 1 to 8 . 

Figures 1 and 2 show the role of magnetic anisotropy for different resonator frequencies. 
Increasing the anisotropy, generally, shifts the delay time, widens the coherence pulse, and 
decreases the coherence peak maximum, and the reversed spin polarization. 

Figures 3 and 4 demonstrate the role of the resonator attenuation. The larger 7 , the 
higher the coherence peak maximum and the shorter the delay time. Spin reversal is better 
pronounced for larger attenuations. 

Figures 5 and 6 illustrate the role of the coupling strength between the magnet and 
resonator. The larger coupling parameter leads to a shorter delay time and higher coherence 
maximum. The value of ca does not influence much the behavior, when the anisotropy is 
weak, A <C 1. Stronger coupling leads to a better spin reversal. 

Figures 7 and 8 show the role of anisotropy for a large value of the coupling parame¬ 
ter. The stronger anisotropy increases the delay time and diminishes the coherence peak 
maximum. Spin reversal is more effective for a weaker anisotropy. 

The general form of the above coupling function looks a bit cumbersome, because of 
which we also check its approximate form, discussed in Section 4, which is given by the 
expression 

a ~ ( 772(1 — Az) (1 — e”'’'*) . 


It turns out that this approximate form is applicable under weak anisotropy, with the 
anisotropy parameters A -C 1, but becomes invalid for larger anisotropies. 


9 Spatial distribution of radiation 

The spatial distribution of radiation for atomic systems can be found in references [22U23] . In 
the case of spin systems, the calculational procedure is similar, with the difference that the 
radiation is produced by moving magnetic moments. The other slight difference is connected 
with the chosen geometry. For atomic systems, one usually considers a cylindric sample, with 
the axis along the 2 ; axis. In the setup, related to a magnetic system, as we study here, the 






magnetic sample is inserted into a coil, with the axis along the x axis, while an external 
magnetic held dehnes the axis, so that the sample axis is orthogonal to the z axis. 

The operator of magnetic moment, related to a j-th spin, can be written in the form 



/ioSj — T k* T koSj , 

(81) 

where 


—. ko ( ■ \ —* 

k 2 7 * 5 k ko^z ■ 

(82) 

The radiation intensity 

in the direction of the vector 

r r 

(83) 


r r 


can be represented |T9] as 

/(n,f) = 2wo7o5^y^y(n)(S'+(t)S'“(t)) , 
o 

where 70 is a natural width, 

7o = ^ ^ /^o^o (^0 = , 

and the system form-factor is 


^Pijin) = — |n X exp{ikon ■ Vij) . 

071 


Here we introduce the unit vector 


jl 1 


and take into account that 




— ^ 1,7 |2 — ,,2 

2 ’ Ih-ol h-o • 


Denoting by d the angle between n and Cz, we have 


|n X = 1 — - sin^ ^ ~ 2 ' 

The form factor (86) becomes 

3 

ifijin) = (1 -h cost'd) exp{ikon ■ Vij) . 

iOTT 

The radiation intensity (84) can be separated into two terms, 

I (n, t) = /inc(n, t) /coh(n, t) , 
the incoherent radiation intensity 

-fine(n,i) = 2^070 Ty(n)(7(*)7(*)) > 


(84) 


(85) 


( 86 ) 
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( 88 ) 


(89) 


(90) 



and the coherent radiation intensity 


J-coh (n, t) = 2a;o7o {t)) , (91) 

where 

(^(n) = cpjj{n) = ^ (1 + cos^^?) . (92) 

Taking into account the identity 

S^S- = S{S + 1) - {S]Y + S] 

and the approximate equality 

which is exact for spin S' = 1/2, as well as for large S —)■ oo, makes it possible to rewrite the 
incoherent radiation intensity as 

-fmc(n, t) = 2 Uq'Jo J2vH{S+{S^}) . ( 93 ) 

j 

Introducing the local functions 

u,(t) ^ 1 {sr(t)} , w,(t) ^ T I {S-(t)) 1 = , 

%(<) = / , ( 94 ) 

and using the semiclassical approximation, we come to the incoherent radiation intensity 

/inc(n,t) = 2a;o7oS'y^9?ij(n)[l + Zj{t)] , (95) 

j 


and the coherent radiation intensity 

/coh(n,t) = 2uo'yoS‘^'^(pij{n)u*{t)uj{t) . (96) 

i¥=j 


Assuming that the radiation wave length is comparable or larger than the linear sample 
sizes, we can resort to the uniform approximation, passing to functions (11), (12), and (13). 
In this procedure, we consider the sum 


i¥=j 


99(n)iV^ 


F{kon) 


1 ■ 

m ’ 


in which 


F{kon) 


1 
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N 



gifcon-rj 


i=i 
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Then we hnd the incoherent radiation intensity 


(97) 


(98) 


line (n, t) = 2uJo'yoNSLp{n) [1 + z{t)] 


(99) 



and the coherent radiation intensity 


/coft(n,t) = 2uJo'yoN‘^S‘^Lp{n)F{kon)w{t) . (100) 

The total radiation intensity, integrated over the spherical angles, is 

I{t) = j /(n,t) dkl{n) = Iinc{t) + Icoh{t) , (101) 

containing the incoherent part 



hncit) = 2uo'yoSN[l + z{t)] 

( 102 ) 

and the coherent part 

Icoh{t) = 2uo'yoS‘^N‘^ipow{t) , 

(103) 

with the shape factor 

990 = / 99(n)F(A;on) dn(n) . 

(104) 

When spins are nniformly distribnted in the sample, fnnction (98) 
an integral over the sample volnme: 

can be represented as 


F{kon) = - / dr 

^ Jv 

The direction vector, in terms of spherical coordinates, is 

n = {sin d cos 99 , sin ■d sin 99 , cosd} . 

Integrating over the sample, we meet the integral 

J sin (^cy/ 0 ? — j cos( 6 a;) dx = ^ ^ 2 '" f 

expressed throngh the Bessel fnnction Ji of the hrst kind. Finally, we obtain 

16 sin^ sin -d cos 99 ) Jlik^Ry/ sin^ -d sin^ Lp + cos^ -d) 


F(A;on) = 


sin^ -d cos^ (p A;oi?^(sin^ -d sin^ (p + cos^ -d) 


(105) 


(106) 


(107) 


Recall that ko = Uq/c = 2n/\. 

Since the magnetic sample is inserted into a coil aligned along the axis x, the radiation 
in the direction is absorbed by the coil. Radiation can be emitted only in the x direction, 
when 

^ ('^ = 5’ '^ = “) ■ 

Then 99 ( 11 ) = 3 /( 167 r), which is only twice smaller than its valne in the 2 : direction. 

If the sample is a very narrow cylinder, or represents a linear chain of spins aligned along 
the axis x, then expression (107) leads to 


4sin^ sin-d cos 99 ) 
sin^ -d cos^ 99 


F{kon) 


(R ^ 0) . 


(109) 



Magnetic dipole radiation nsnally corresponds to long waves, snch that the wavelength 
A is longer that the sample linear sizes. To consider this limit, we notice that the Bessel 
fnnction J^^x), at small a; <C 1 , has the asymptotic property 

^ , 1 /X\'^ 1 /X \^'+2 

^ ~ r(i/ + 1 ) V 2 y “ r(i/ + 2 ) ( 2 / 

becanse of which 

X 

M^) - 2 (a: -)■ 0 ) . 

Hence, in the long-wave limit, 

,, , / 27rR ttL a , , 

F{kon) ~ 1 ( ^ ^ 

In this way, althongh a part of radiation is absorbed by the coil, nevertheless, an essential 
part of radiation can be emitted throngh the sides of the sample along the coil axis x. 


10 Discussion 

We have considered a large class of magnetic materials, composed of particles interacting 
throngh exchange interactions. Snch materials are well represented by an anisotropic Heisen¬ 
berg Hamiltonian. By conpling a magnetic sample to a resonant electric circnit it is possible 
to efficiently regnlate spin dynamics. The present investigation complements the earlier 
stndies of spin dynamics in different materials, accomplished for polarized nnclei, magnetic 
nanomolecnles, and magnetic nanoclnsters. The systems composed of these particles pos¬ 
sess several common properties. Snch systems, except polarized nnclei, exhibit magnetic 
anisotropy. For example, the spins of magnetic nanomolecnles are frozen below the block¬ 
ing temperatnre Tb ~ 1 — 10 K, with the frozen magnetization protected by an anisotropy 
barrier of energy Ea ~ 10 — 100 K. For magnetic nanoclnsters, the blocking temperatnre 
is Tb ~ 10 — 100 K. Polarized nnclei can be represented by polarized protons in hydro¬ 
genated materials, snch as propanediol, bntanol, and ammonia. Althongh these materials 
do not enjoy magnetic anisotropy, however their spins, at low temperatnre, can remain po¬ 
larized for extremely long time. The common featnre of all above materials is that their spin 
interactions are described by dipolar forces. 

However, there exists a very wide class of materials made of particles with magnetic 
exchange interactions. Also, recently there have appeared a new type of magnetic nanoma¬ 
terials characterized by exchange interactions, snch as magnetic graphene, where magnetic 
properties are indnced by defects [H]. 

Graphene is a two-dimensional carbon material intermediate between an insnlator and a 
metal pn[25] . Carbon-carbon spacing is a ~ 1.42 A, and its snrface density is p ~ 3.9 x 10 ^® 

_9 

cm . 

Magnetic graphene can be well described by an anisotropic Heisenberg model [2H1]. For 
defects on a zigzag edge, one has an exchange interaction potential J ~ 0.1 eV, or J ~ 10“^^ 
erg. This gives J/h ^ 10^^ s“^ and J/ks ~ 10^ K. Dipolar interactions are mnch weaker, with 
pl/a^ 10 srg, or jxj^/[h(E) ~ 10^*^ s In temperatnre nnits, this gives [kBd^) ~ 
0.1 K. Magnetic anisotropy is not strong, with AJ/J ~ 10“^, hence, AJ/h ~ 10^° s“^. 
Thns, 72 ~ ^%/{ha^) ~ 10^° s“h If 73 ~ AJ/h, then 73 ~ 10^° s“^, that is, of the same 





order as 72 . In the case of an external magnetic field -Bq = 1 T, the Zeeman frequency is 
<^o = Ih-B-Bol/^ 10^^ s Therefore the radiation wavelength is A ~ 10 cm. Since there 
are many ways of generating defects in graphene, the system parameters can be varied. 

It is important to stress that a self-organized spin dynamics, when there is no initial 
coherence imposed onto the sample, cannot be correctly described by phenomenological 
equations, like Landau-Lifshitz, Gilbert, or Bloch equations. This is why, we use a micro¬ 
scopic approach with a realistic Hamiltonian. Such an approach makes it possible to take 
into account quantum fluctuations, crucially important at the beginning of spin relaxation. 
Also, the use of a microscopic model has an advantage of containing well dehned parameters 
associated with the considered Hamiltonian. 

When the system is prepared in a nonequilibrium state, and no initial coherence is im¬ 
posed onto the sample, spin wave fluctuations serve as a triggering mechanism starting spin 
motion. It is possible to show that spin wave fluctuations is the sole triggering mechanism, 
while thermal Nyquist noise of the coil cannot play the role of a trigger. 

In the process of motion, spins are collectivized by means of the resonator feedback 
held. Superradiance in spin systems cannot be caused by photon exchange. That is, spin 
superradiance is due to the Purcell effect and is impossible without a resonator. This is 
contrary to atomic systems, where superradiance develops as a Dicke effect. 

The microscopic equations of motion are investigated by employing stochastic mean- 
held approximation and scale separation approach. The equations for guiding centers are 
solved numerically. Depending on initial conditions and system parameters, there can exist 
diherent regimes of spin dynamics, slow free relaxation during the time Ti, free induction in 
time T 2 , weak superradiance in time slightly shorter tan T 2 , pure superradiance and triggered 
superradiance, when the reversal time of magnetization can be made much shorter than T 2 , of 
order 10“^^ s or 10“^^ s. The regime of punctuated superradiance can be realized, producing 
a sequence of coherent pulses. 

We have analyzed the spatial distribution of radiation produced by moving magnetic 
moments. Although a part of radiation is absorbed by the coil surrounding the sample, 
anyway, radiation can be emitted through the sides of the sample, where there is no coil. 

The studied effects can be used in a variety of problems in spintronics and in quantum 
information processing. 
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Figure Captions 


Figure 1. Role of magnetic anisotropy in spin dynamics for the parameters 7 = 1 , 
g = 10, and u = 10, under different anisotropy parameters, R = 0.1 (solid line), A = 0.5 
(dashed line), and A = 1 (dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 


Figure 2. Role of magnetic anisotropy in spin dynamics for the parameters 7 = 1 , 
g = 10, but u = 100, under different anisotropy parameters, R = 0.1 (solid line), A = 0.5 
(dashed line), and A = 1 (dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 


Figure 3. Role of resonator attenuation for the parameters 5 ^ = 10, R = 0.1, and oj = 10, 
under different attenuation parameters, 7 = 1 (solid line), 7 = 10 (dashed line), and 7 = 100 
(dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 


Figure 4. Role of resonator attenuation for the parameters g = 10, A = 0.1, but 
u = 100 , under different attenuation parameters, 7 = 1 (solid line), 7 = 10 (dashed line), 
and 7 = 100 (dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 


Figure 5. Role of magnet-resonator coupling for the parameters 7 = 10, A = 0.1, and 
u = 10 , under different coupling parameters g = 10 (solid line) and g = 100 (dashed line): 
(a) Coherence intensity; (b) Spin polarization. 


Figure 6 . Role of magnet-resonator coupling for the parameters 7 = 10, R = 0.1, but 
u = 100 , under different coupling parameters, g = 10 (solid line) and g = 100 (dashed line): 
(a) Coherence intensity; (b) Spin polarization. 


Figure 7. Role of magnetic anisotropy for the parameters 7 = 10, 5 ^ = 100, and 00 = 10, 
under different anisotropy parameters, A = 0 (solid line), R = 0.1 (dashed line), A = 0.5 
(dashed-dotted line), and A = 1 (dashed line with dots): (a) Coherence intensity; (b) Spin 
polarization. 


Figure 8 . Role of magnetic anisotropy for the parameters 7 = 10, 5 ^ = 100, but u = 100, 
under different anisotropy parameters, A = 0 (solid line), A = 0.1 ( dashed line), A = 0.5 
(dashed-dotted line), and A = 1 (dotted line): (a) Coherence intensity; (b) Spin polarization. 





Figure 1: Role of magnetic anisotropy in spin dynamics for the parameters 7 = 1, 5 ^ = 10, 
and oj = 10, under different anisotropy parameters, R = 0.1 (solid line), A = 0.5 (dashed 
line), and A = 1 (dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 




Figure 2: Role of magnetic anisotropy in spin dynamics for the parameters 7 = 1, 5 ^ = 10, 
but OJ = 100, under different anisotropy parameters, A = 0.1 (solid line), A = 0.5 (dashed 
line), and A = 1 (dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 














Figure 3: Role of resonator attenuation for the parameters g = A = 0.1, and oj = 10, 
under different attenuation parameters, 7 = 1 (solid line), 7 = 10 (dashed line), and 7 = 100 
(dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 




Figure 4: Role of resonator attenuation for the parameters g = IQ, A = 0.1, but u = 100, 
under different attenuation parameters, 7 = 1 (solid line), 7 = 10 (dashed line), and 7 = 100 
(dashed-dotted line): (a) Coherence intensity; (b) Spin polarization. 
















Figure 5: Role of magnet-resonator coupling for the parameters 7 = 10, A = 0.1, and 
u = 10 , under different coupling parameters g = IQ (solid line) and g = 100 (dashed line): 
(a) Coherence intensity; (b) Spin polarization. 




Figure 6 : Role of magnet-resonator coupling for the parameters 7 = 10, A = 0.1, but 
oj = 100 , under different coupling parameters, g = IQ (solid line) and g = 100 (dashed line): 
(a) Coherence intensity; (b) Spin polarization. 






















Figure 7: Role of magnetic anisotropy for the parameters 7 = 10, 5 ^ = 100, and u = 10, 
under different anisotropy parameters, A = 0 (solid line), R = 0.1 (dashed line), A = 0.5 
(dashed-dotted line), and A = 1 (dashed line with dots): (a) Coherence intensity; (b) Spin 
polarization. 




Figure 8 : Role of magnetic anisotropy for the parameters 7 = 10, = 100, but oj = 100, 

under different anisotropy parameters, ^4 = 0 (solid line), A = 0.1 ( dashed line), A = 0.5 
(dashed-dotted line), and A = 1 (dotted line): (a) Coherence intensity; (b) Spin polarization. 


























